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pi 1 The paper develop a new approach to the justification of Gibbs canonical 

. distribution for Hamiltonian systems with finite number of degrees of freedom. 

It uses the condition of nonintegrability of the ensemble of weak interacting 
Hamiltonian systems. 
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1. Gibbs distribution. We consider the probability distribution in the phase 
space of Hamiltonian system with the density 

p = ce" 7 ^, (1) 



where if is a Hamiltonian, r is an absolute temperature, k is the Boltzmann con- 



stant, c is a normalized factor. It plays the key role in the equilibrium statistical 
mechanics. Gibbs show in [1] that the averaging with respect to probability measure 
with density (1) give rise to the fundamental relations of equilibrium thermodynam- 
ic . ics. 

To deduce the canonical Gibbs distribution one usually consider the ensemble of 
Hamiltonian systems with Hamiltonian function of the following form 



> 

•i-H . 

^ ■ where 



H = H (P,Q)+eH 1 (P,Q), (2) 



s=l 



p« = (pi' ) ) ...,pW) ) q ^ = (q[ s \...,q£). (3) 

Thus at s = the system with Hamiltonian (2) is decomposed on n identical systems 
with m degrees of freedom and Hamiltonian H . The canonical variables P, Q are 
the momenta p^ 1 ', . . . , and coordinates g^ 1 -*, . . . , of separate subsystems. 
The perturbing function Hi is the energy of interaction of n subsystems; it usually 
depends on their coordinates Q. Small parameter e is the characteristic of intensity 
of subsystems' interaction. 
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We consider the case, when the Hamiltonian 7i is sufficiently smooth with re- 
spect to variables P, Q. However, in application we often see cases with singular 
interaction. The classical example is the Boltzmann-Gibbs gas, the assembly of rigid 
balls in cube that elastically collide with each other (see [1,2,3]). 

The traditional approach to the deduction of Gibbs distribution suggested by 
Fowler and Darwin ([4], the rigorous exposition see in [5,6]) essentially uses the 
ergodic hypothesis: for all small e > the Hamiltonian system with Hamiltonian (2) 
is ergodic on fixed energy manifolds 7i = const. With some additional conditions 
some systems with Hamiltonian Hq are distributed in accordance with formula (1) 
as e — > and n — > oo. 

But the proof of ergodic hypothesis for specific Hamiltonian systems is usually 
pretty difficult problem. Moreover the ergodic hypothesis often contradicts with 
results of KAM theory. In particular if the Hamiltonian system with Hamiltonian Hq 
is completely integrable and the energy surfaces Hq = const are compact, then the 
ergodic property is not present with certainty. 

In view of this remark we can set the following interesting problem: prove that 
in analytical (or even in infinitely differentiable) case if the energy surfaces H = 
const are compact, then the Hamiltonian system with Hamiltonian (2) never satisfy 
ergodic hypothesis. 

We can try to modify the Fowler-Darwin method assuming that e ^ and n 
depend on e in such way that n(e) — > oo as e — > 0. We can assume that for some 
functions e i— > n(e) the system with Hamiltonian (2) is ergodic on the surfaces 7i = 
const as a result of huge number of its degrees of freedom. This somewhat weakened 
version of ergodic hypothesis is closely related to the unsolved problem of estimation 
of small parameter in KAM theory, when the "last" Kolmogorov torus disappears. 
The weaker conjecture on transitivity of system with Hamiltonian (2) on energy 
surfaces 7i = const for large values of n and small e is not proved yet. If this 
problem has the positive answer, then we can assert the presence of diffusion in 
Hamiltonian systems with many degrees of freedom (see [7-9]). 

2. Probability density as an integral of Hamiltonian equations. A different 
approach to the deduction of canonical Gibbs distribution was proposed in the paper 
[10]. This approach is based on the fact that the stationary density of probability 
distribution is the integral of Hamiltonian differential equations uniquely defined in 
the whole phase space. In this case the number of interacting subsystem n > 2 is 
fixed. 

More exactly in [10] we consider the case, when the subsystems have one degree 
of freedom: m — 1. Under some natural condition we can proceed to the angle- 
action variables in each subsystem, and using the well-known Poincare method we 
can obtain the constructive condition of nonexistence of new integrals (see [11, 12]). 
The results of paper [10] can be easily converted to the more general case, when Hq 



2 



is a Hamiltonian of completely integrable system. We are going to find out the 
sufficient conditions (constructive if possible) of nonintegrability of systems (2)-(3). 

We study the conditions of existence of an integral T(P, Q, e) of the canonical 
differential equations 

with a Hamiltonian 7i of the form (2)-(3). We emphasize that the integral JF 
depends on the parameter e. Poincare considered the analytical case; in particular 
we can construct the integral T as a power series with respect to e. We suppose 
that T is a function of class C 2 with respect to all the variables P, Q and e. Hence 
we can suppose 

F = F (P,Q)+eF 1 (P,Q) + o(e), (5) 

where JF and T\ are functions of class C 2 and C 1 with respect to P and Q corre- 
spondingly. Probably we can weaken the requirements to the class of smoothness of 
integral T, and the following arguments still will be correct. But this is a separate 
problem and we are not going to discuss it here. 



3. Non-perturbed problem. Suppose e — 0. Then we have a system of n 
independent subsystems. It is strongly nonergodic: at e = system of differential 
equations (4) has n independent first integrals 

H s = H (pM, q {s) ), 1 < s < n. (6) 

It is clear that the function JF from expansion (5) is a first integral of this unjointed 
system. Let's show that under the specific conditions the function JF depends 
only on Hi, . . . , H n . In particular these conditions will imply that any separate 
subsystem with m degrees of freedom does not nave first integrals independent on 
the integral of energy. The ideas of our arguments follows Poincare method [11]. 

Let M be a phase space of Hamiltonian system with Hamiltonian function (6). 
Certainly for all s these spaces are identical. A phase space of new system is the 
direct product 

M = M x . . . x M, dim.M = 2nm. 
Let h s be a value of total energy of system with number s, and 

J>.) = {p (s) , q {s) : H s (p^, ,«) = h s } 

is the a energy surface. If the value of h s is uncritical, then is a smooth 2m — 1- 
dimensional manifold. At fixed values of h — (hi, ... , h n ) and e = unjoined 
Hamiltonian system (4) is reduced to the direct product of dynamical systems de- 
fined on 

S(h) = J2(hi) x . . . xJ2(h n ) C M. 
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The ergodic property of system with Hamiltonian H s on ^2(h s ) does not neces- 
sarily imply the constancy of integral JF on the invariant set S(h). Consider the 
following simple example 

EXAMPLE. Suppose the following dynamical system 

±i = tvi, yj = ujj; i, j — 1, k (7) 

with constant incommensurable systems uj = (u>i, . . . , u>k) is defined on direct prod- 
uct of /c-dimensional tori T h {x mod 2n} x T k {y mod 2tt}. According to the Weyl 
theorem, each separate subsystem is ergodic on T h . But equations (7) have single- 
valued nonconstant integrals sin(a;i — Hi) (1 < % < k). 

Remark. However, if Hamiltonian systems are weakly mixing (mixing) sys- 
tems on ^2(h s ), then their direct product also has weakly mixing (mixing) property 
on S(h). In particular in these cases the function JF is constant on any connected 
component of manifold S(h). 

Let T 1 be a nondegenerate periodic trajectory of system with number s with 
energy h s , T s its period, and uo s = 2n/T s its frequency According to Floquet- 
Lyapunov theorem, in a neighborhood of this trajectory on ^2(h s ) we can express 

(s) (s) 

the coordinates (p s mod 2ir, z\ , . . . , z\^_ 2 ^ so that in the new variables the motion 
equation obtain the following form: 

up s = cu s + f 8 (<p 9 , z^), z^ = Q s z^ + g s (p s , z^). (8) 

Here f s = 0(\z^\), g s = o(\z^\), and constant square matrix il s of order 2m — 2 
is nondegenerate. Assuming in (8) z^ = we obtain the equation of periodic 
trajectory: 

Ps=u s (1 < s < n). (9) 

According to the assumption on non degeneracy of periodic trajectory T 1 , non- 
degenerate periodic trajectories with similar period are situated on close energy 
surfaces ^2(h s ); periods and frequencies continuously depend on energy h s . 

It is clear that the direct product T 1 x . . . x T 1 = T™ is n— dimensional invariant 
torus of canonical system of differential equations (4) at e = situated on S(h). 
In the neighborhood of this torus the equation of motion have form (8). Hence, 
such torus is reducible and nondegenerate (see, for example, [12]). On the torus the 
equation are reduced to a conditionally-periodic form (9). As usually, we call an 
invariant torus nonresonance if the frequencies oui , . . . , ou n are independent on the 
ring of integers. In the further analysis the following condition is essential 

A) For almost all admissible values of h G M ra nonresonance tori are everywhere 
dense on the manifold S(h). 

EXAMPLE. Let separate subsystems describe the inertial motion on the man- 
ifold Af of negative curvature. Energy h is non-negative. All periodic trajectories 
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with positive energy are hyperbolical; hence they are not degenerate. Periodic tra- 
jectories are the motions on closed geodetics on M with different speed. If / is a 
length of closed geodetic, then the period is equal to Hence the frequency cu is 
defined by the formula /ttt 



Since the lengths of n geodetics l\, . . . , l n are fixed, then for almost all positive 
values of energy hi, . . . , h n the frequencies u>i, . . . , uj n are incommensurable. It is 
possible to show (and this is a separate problem) that in a considered situation 
condition A is fulfilled. 

Proposition 1. // condition A is fulfilled, then for all h G M n the function JF 
is constant on any connected component of S(h). 

Proof. Since JF is an integral of system of equations (4) at e — 0, then (according 
to the Kronecker theorem) JF is constant on any nonresonance torus T n . Since this 
torus is reducible and nondegenerate, then dT§ = in points T (see [12], ch. IV). 
According to condition A, for almost all values of h G MP nonresonance tori are 
everywhere dense on S(h). Hence, dJ-'o = on such manifolds S(h). Therefore T§ 
is constant on their connected components. For other values of h the conclusion of 
proposition 1 follows by continuity. 

Remark. The proof shows that in condition A instead of "for almost all admis- 
sible values of h G M. n " we can say "for everywhere dense set of values of h G M ra ". 
However, such weakening of condition A practically does not give anything new. 

In further analysis we will use the proposition on everywhere density of the set 
of maximum resonance tori (when all frequencies are rationally expressed through 
one frequency). This condition together with condition A produces "an alternation" 
of resonance and nonresonance invariant tori and replaces the condition of nonde- 
generacy of non-perturbed completely integrable system in the Poincare theory. 

4. Energy surfaces. Let us consider a case, when function the H : M — > M. 
has a finite number of critical values a\ < a 2 < . . . < a r , and ai = miniJ . Such 
situation is common in applications. When the total energy h passes through the 
critical value, the continuous dependence of energy surface XX^o) on ho is lost. In 
that moment its topology generally changes. 

In fig. 1 we present the plot of Hamiltonian H Q with three critical values. The 
points a 1 and a 3 are stationary values of H , and the critical point a 2 is not a sta- 
tionary value. The presence of the nonstationary critical points is the characteristic 
property of potentials describing gravitational or Coulomb interaction. 

Let's denote as K ili2 ___ in an open parallelepiped in M, n = {hi, h n }. This 
parallelepiped is a direct product of the intervals 




a«i < hi < a il+ i, . . . , a,i n < h n < a in+ i. 



(10) 
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If number i s + 1 is larger than r, then we replace Oj s+ i with a symbol oo. In fig. 2 
these domains are shown for n = 2 and r = 3. Each point h G -Kj, % = (ii, . . . , i n ) 
corresponds to a smooth regular manifold, which may consist of several connected 
parts. The quantity of connected components of S(h) does not depend on a point 
h G Ki] we denote this number by symbol Xj. 

Let us introduce in the phase space M. open areas defined by the inequalities 
similar to (10): 

a h < < a il+1 , . . . , a in < H n {p (n \ q {n) ) < a in+1 . 

It is clear that the closure of these domains in the whole covers all M.. Also each 
/Q has exactly Xj connected components. 

Proposition 2. For any connected component of domain /Cj there exists the 
continuously differ entiable function 

such that the following equality is fulfilled in this domain 

Fo = f i (H 1 ,H 2 ,...,H n ). (11) 

Remark. Actually function fa belongs to the class of smoothness C 2 . However 
it is not essential for the further analysis. 

Proof. It is clear that the domain /Q is foliated to the regular surfaces S. The 
function JF is constant on these surfaces (more exactly on their connected com- 
ponents) (the proposition 1). Hence, on any connected component /Q the func- 
tion To has natural representation (11). By definition in any point /Q functions 
Hi, . . . , H n are independent. Therefore we can introduce locally new coordinates 
so that Hi, . . . , H n will appear as n of new variables. The transition to such co- 
ordinates is carried out with certainty with the help of continuously differentiable 
reversible transformation. In new variables the function JF is continuously differen- 
tiable and depends only on n variables Hi, . . . , H n . The proposition is proved. 

5. Resonances. Let us consider again an invariant torus T of non-perturbed 
system. The equations of motion on the torus are reduced to form (9). This torus 
we call completely resonance if there exist n — 1 linearly independent integer vectors 

U= (Ui, U n ), V = (Vi, V n ), W = (Wi, . . . , w n ), 

such that 

(u, uj) = UiCUi + . . . + u n oj n = 0, (u, uj) — 0, . . . , (w, uj) = 0. (12) 
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In other words all frequencies u s are rationally expressed through one of them. This 
is equivalent to the proposition that all solutions of differential equations (9) on the 
torus T™ are periodic with the same period. 

Let $ be a restriction of perturbing function 7i : M — » R on the invariant torus 
T n . It is clear that $ is 2n— periodic function of (</?i, . . . , ip n ) = (p. We associate it 
with multiple Fourier series 

$ = ^ (p k expi(k, <p). (13) 

A nonresonance torus we call the Poincare torus if the factors of Fourier decom- 
position (13) with numbers u, v , . . . , w are nonzero. Since the Poincare tori consist 
of the separate closed trajectories, they have no "rigidity" property and collapse af- 
ter the addition of perturbation. We can not exclude the possibility that the family 
of degenerate periodic trajectories, components of the Poincare torus, at perturba- 
tion give rise to the finite number of nondegenerate periodic trajectories with close 
transition. 

Let us introduce finally the Poincare set P C R n . It is a set of points from 
R n = {hi, . . . , h n }, which are the images of the Poincare tori under the "energy" 
mapping M. — > R n : a point with coordinates 

(P,Q) = (p (1) ,...,p (b) , q (n \ ...,g (n) ) 

passes to a point with coordinates 

hi = Hi{pM,qV), K = H n tf"\qV>). 

Remark. In the Poincare theory [11, 12] the usually supposition is that the non- 
perturbed system with the Hamiltonian TLq is completely integrable and nondegener- 
ate. Therefore the set of its first integrals is a set of action variables "numerating" 
the invariant tori. The Poincare set (introduced in [12]) is defined here as a set 
of points in space of action variables corresponding to the completely resonance tori 
collapsing after the addition of perturbation. In our case functions Hi, . . . , H n are 
the set of integrals of the non-perturbed problem and the Poincare set is a set of 
points in space of values of these integrals. 

Proposition 3. In points of the Poincare set P functions 

n 

n o = ^2 and = f(H u . . . , H n ) 

8 = 1 

are dependent. 

Proof. Let { } be a Poisson bracket connected with symplectic structure on M.. 
Since function T is the first integral of initial system (4), then {H, J 7 } = for all 
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values of e. Since {H , JF } = 0, then 

lim {H ' ^ = 0. (14) 
Using decomposition (2) and (5) we receive from (14) the equality 

{Ho, ^i} = {-Fo, Wi}. (15) 

Well known that {JF , H\\ is a derivative of T\ by virtue of the initial system 
of differential equations with Hamiltonian Ho. The Poisson bracket {JF , Hi} has 
similar sense. Besides by formula (11) we obtain 

{^o, Wi} = £ TS^-W. Wi>- ( 16 ) 

8 = 1 S 

Now let us restrict equality (15) on the invariant torus T n . It is clear that 
{Ho, Fi} is equal to the derivative of restriction of function T\ on T n by virtue of 
system of differential equations (9). Let \I> : T" — > R be a restriction of JF X on T n 
and 

* = ^expi(k, <p) (17) 

its Fourier series. In points tp G Z n the left part of (5.4) becomes 

j(fc, u)tp k exp j(fc, (/?). (18) 



In view of formulas (13) and (16) the right part of relation (15) in points of the 
invariant torus is equal to 

^ i \^Jf klUJl + ■■■ + -M-k n ^nj (fkexpi(k, <p). (19) 

Comparing (18) and (19) we obtain a chain of equalities 

(k, u)^ k = (f fciWi + . . . + £ n k n uJn)ipk, k G Z n . (20) 

Now let T n be the Poincare torus. Setting k equal to w, v , ...,wwe obtain that 
w as a vector of n— dimensional space is orthogonal to the hyperplane II generated 
by linearly independent vectors u, v, . . . , w. Since (u, u) — . . . — (w, u>) = 0, and 
ip u ^ 0, . . . , <f w ^ 0, then we obtain from (20) n — 1 linear relations: 

ui(fiWi) + • • • + Un(fnW„) = 0, 

tWl(Cl^l) + • • • + Wn(^n^n) = 0. 
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Hence the vector with components £iUi, . . . , £ n cj„ is orthogonal to the hyperplane II 
and consequently it is collinear to the vector u = (ui, . . . , ou n ). Hence £1 = . . . = £ n . 

Thus in points of the Poincare set the derivatives . . . , -J^f - are equal to each 
other. It means obviously the dependence of functions Ho and To- 

Remark. Since functions $ and ^ are assumed to be only once continuously 
differentiable, then the Fourier series (13), (17)-(19) may diverge. In this case 
relations (20) can be deduced in another way. For this purpose we multiply the 
derivative with respect to time 

■ dS d<h 

$ = JT~ + • • • + 7T^™ 
dipt dip n 

on expi(k, ip), apply the operation of averaging overT n 

j^Y J (-)^i 

and integrate by parts. In result we obtain the left part of relation (20) up to the 
factor —i. The right part of (20) is obtained similarly. 
Let us introduce the additional condition 

B) The Poincare set P is everywhere dense in any parallelepiped Ki. 
The proposition 3 imply the following corollary 

Corollary. // conditions A and B are fulfilled, then the functions To and JF are 
everywhere dependent. 

Hence we immediately obtain the following representation on any parallelepiped 

To = F t (Ho), (21) 

where Fj is some continuously differentiable function. 

Indeed, introducing the new variables Hi, H n -i,TLo = H s instead of 
Hi, . . . , H n we can write down formula (11) in the another form: 

— fi(H±, . . . , H n _i, Ho — Hi — ... — H n _i). 

Since functions To and Ho are dependent, then the right part of this equality actually 
does not depend on Hi, . . . , H n _i. 



6. Deduction of Gibbs distribution. Now let e tend to zero. At the limit we 

obtain the n unjoined subsystems moving independently: the change of the initial 
data p( l \ for I ^ s does not affect the dynamics of subsystem with number s. In 
order to be consistent we shall also assume at e = that the subsystem with number 
s being m some fixed state (p (s) , g (s) ) € M is a random event. The following natural 
condition plays the main role in the process of deduction of Gibbs distribution 
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C) These random events are independent. 

n is a density of probability distribution of the subsystem 

with number s, and 

p(P,Q, e) = p (P,Q)+O(e) 

is a density of probability distribution in the initial system with Hamiltonian (2), 
then using condition C and the rules of multiplication of probabilities of independent 
events we obtain as e — > the following equality: 

Po = pi ... p n - (22) 

Equality (22) is also called the Gibbs hypothesis on the preservation of thermo- 
dynamic equilibrium of subsystems at vanishing interaction ([1], see also [13]). The 
sense of this term will be explained below. 

Our main result is the following theorem 

Theorem. Suppose conditions A, B and C are fulfilled. Then 

p s = c s e - "^", c s = const > (23) 

for all 1 < s < n. 

In particular, according to (22), 

Po = c e~"^, c = c 1 . . . c n . 

From (22) we see that all separate subsystems have the same distribution. We 
compute the factors c s using the normalizing condition 

j Ps d n pd n q=l. 

M 

Proof of theorem. First note that p s is a function of Hamiltonian H s only, and 
it is continuously differentiate in all open intervals 

(ai, a 2 ), (a 2 , a 3 ), (a r , oo). (24) 

More exactly the number of such functions is equal to the number of the connected 
components of level surface ^2(h s ), when the energy h s changes in each of intervals 
(24). Some of these functions may coincide. 

Indeed, for almost all h s > a\ the Hamiltonian system with Hamiltonian H s 
has everywhere dense set of nondegenerate periodic trajectories on energy surfaces 
^2(h s ). Otherwise condition A is not fulfilled because of the identity of separate 
subsystems. Then, according to Poincare [11], the points of nondegenerate periodic 
trajectories are stationary for the restriction of any first integral on ^2(h s ). The 
continuity imply that the first integrals of Hamiltonian system with Hamiltonian 
function H s are constant on the connected components of ^(h s ). At last we should 
note that p s is the first integral and use the (simplified) arguments of section 4. 
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Now let us consider again equation (22) true on any parallelepiped Kf 

p (#i + ... + H n )= p^Hi) . . . Pn (H n ). (25) 

This functional equation is easily solved. We differentiate (25) sequentially with 
respect to Hi, . . . , H n and divide the result on the product p± . . . p n . In result we 
obtain the following chain of equations 

P± = ... = PjL = -ft. 
Pi Pn 

Here f3 is some constant independent on the number s. Hence 

p s = c s e~ ~P Ha , c s = const. (26) 

The dimension of constant (5 is equal to the inverse energy dimension. Usually 
one suppose that f3 = (Arr)" 1 , where r is the absolute temperature, and k is the 
Boltzmann constant. 

We should note that formula (26) may depend on the choice of multiindex % = 
(ii, . . . , i n ). More precisely, any connected component of the set /Q C M. has its 
own set of the factors (3 and c s in (26). 

However, we can easily show that the constant (3 has an universal character. 
Indeed, suppose the constants j3 in formula (26) are equal to the values f3\, . . . , f3 n 
on some connected component of domain /Q with some index %. Then, according to 
(6.1), in this domain 

p = Coe £/^. (27) 

If some P s are not equal to each other, then functions (27) and 7i = J2(H S ) are 
independent. But this statement contradicts to the corollary of proposition 3. 

This argument has the evident physical meaning: at thermodynamic equilibrium 
all components of system have identical temperature. 

The last remaining possibility is that constants c s in formula (26) are different 
on the different intervals of values of energy (24). But in reality this possibility is 
not realized because of the continuity property of functions p s : M — > R. 

The theorem is completely proved. 

In conclusion of the work we shall make one important remark. The Gibbs 
theory presented in [1] does not imply that the densities of probability distributions 
pi, . . . , p n should be continuous functions on M. We can consider more general 
situation and assume, for example, that functions p s are continuously differentiate 
only on those domains of phase space M, in which energy is contained between its 
neighboring critical values 

a r < H s < a r+ i (r = 1, . . . , p; a p+1 = oo). (28) 
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Naturalness of such assumption is already evident if we consider the example of 
mathematical pendulum: the separatrices on phase cylinder separate the domains 
with essentially different type of motion (fluctuations and rotations). 

Applying the developed above method we again obtain formula (23), but the 
constants c s have different values in different domains (28). Moreover, these con- 
stants may be different for different connected components of domains (28). It is 
quite possible that such generalized discontinuous Gibbs distribution could be useful 
for the study of concrete thermodynamic systems. 

Let us assume, for example, that the phase space M has only two domains M + 
and M_ of form (28). In domains M± we have the following densities of distributions 

H 

c±e kT . 

We calculate the constants c + and c_ using the normalizing condition 

c + J e~£d n pd n q + c- J e~£ d n pd n q = 1. 

M + M- 

If the difference A = c + — c_ is given, then factors c± are uniquely defined by this 
equality: 

c+ J e~£d n pd n q = 1 + A J e-^d n pd n q, 

M M- 

c J e-£d n pd n q=l-A J e~£d n pd n q. 

M M+ 

Since c± > 0, then the right parts of these equalities are positive. It happens with 
certainty if the jump A satisfies the following inequality 



IAI < 



e krd n pd n q 



M 



-1 



Could we obtain the probability distribution with piecewise smooth function of 
distribution using the Fowler-Darwin method? 

The paper is prepared with the financial support of RFBR(99-01-0196) and IN- 
TAS. 
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